Cavity solitons in bidirectional lasers 
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We show theoretically that a broad area bidirectional laser with slightly different cavity losses for 
the two counterpropagating fields sustains cavity solitons (CSs). These structures are complemen- 
tary, i.e., there is a bright (dark) CS in the field with more (less) losses. Interestingly, the CSs can 
be written/erased by injecting suitable pulses in any of the two counterpropagating fields. 



Cavity solitons are localized structures that form in 
broad area nonlinear optical cavities in the plane orthog- 
onal to the field propagation direction. They are inter- 
esting not only from a fundamental point of view but also 
because their potentialities in optical information storage 
and processing [lj]. 

Interestingly, the most important nonlinear optical 
cavity, namely the laser, does not exhibit CSs easily, i.e., 
something must be done for usual lasers display localized 
structures different from light vortices. One possibility 
is introducing a suitable saturable absorber in the laser 
cavity 0. Another possibility consists in rocking the 
laser Q, i.e., in injecting an amplitude modulated signal. 
These two ways have been experimentally demonstrated 
in lasers Q or laser-like systems @ . Apart from these, it 
has been numerically shown that cascade lasers [||, and 
lasers with a dense active medium Q , can sustain CSs (in 
the second case some intracavity spatial filtering must be 
done), but up to now there is no experimental evidence 
for these predictions. 

In this letter we numerically show that CSs can form 
in bidirectional lasers provided that the two counterprop- 
agating fields experience slightly different losses (up to a 
few per cent). Interestingly, CSs in the two fields are com- 
plementary (a bright CS in one field complements with 
a dark CS in the other field) and can be written/erased 
by acting on only one of the two fields. 

Consider a broad-area ring-cavity filled with a homo- 
geneously broadened two-level active medium which is 
incoherently pumped. The Maxwell-Bloch equations of 
this system, within the uniform field limit and rotating 
wave-approximations, are Q 

d t E n = -K n E n -iA (P n ) +iaV 2 L E n , (la) 

d t P n = -7_l(1-*A)P„ (lb) 

+i 1± D {E n + E m exp [(-1)" 2ikz]} 

d t D = - 7| | [D - 1 - Im (Ej.Pt + E 2 P;)] , (lc) 

with n, m — 1,2 and n =/= to. In the above equa- 
tions E n (x,y,t), denote the field amplitudes of the two 
(n = 1, 2) counterpropagating fields traveling along the z 
direction, P n (x,y,z,t) are their respective medium po- 
larization sources, D (x, y, z, t) denotes the population in- 



version, and 

(P n (x,y,t)) = -|- J dzP n (x,y,z,t) , (2) 
o 

are the axially-averaged medium polarizations with L z 
the laser cavity length. = (9^ + dy) is the trans- 

verse Laplacian and regarding the parameters, A = 
(luc — ^at) /7j_ is the normalized cavity detuning (loc 
and oj at are, respectively, the cavity and atomic tran- 
sition frequencies), A is the pump parameter, a is the 
diffraction coefficient, and the field, medium polarization 
and inversion decay rates are denoted by n n , and 711, 
respectively. 

In order to keep the mathematical modelling as simple 
as possible, we shall restrict our study here to class-A 
lasers (i.e., those verifying 7^ 3> 711 ki,k 2 ) and shall 
consider perfect resonance, i.e., A = 0. For class-A lasers 
all material variables can be adiabatically eliminated. Af- 
ter taking d t P n — d t D — 0, one easily gets 

i{E n + E m exp [(-!)" 2ikz}} 
1 + |£i | 2 + \E 2 \ + E x E* 2 e^ kz + c.c. 

Now we assume that the pump parameter value is close to 
its value at lasing threshold, so that field amplitudes are 
small quantities. Then, after Taylor expanding Eq. ([3]), 

it follows that (P n ) = iE n (l - \E n \ 2 - 2 \E m \ 2 ^ , and 

then Eqs. (p} reduce to 

d T F 1 = (A-l- |F 1 | 2 -2|F 2 | 2 )F 1 +iVl^ 1 ,(4a) 

d T F 2 - (A -a- \F 2 \ 2 -2\F 1 \ 2 )F 2 + iVlF 2 {Ah) 

where we have introduced F„ = V AE n , Vj^ = aV^, 
t = Kit, and o~ — k 2 /ki being £ = xj^fa the dimen- 
sionless transverse coordinate. This is the model for a 
resonant class-A bidirectional laser we shall study here, 
which consists of a pair of coupled Ginzburg-Landau 
equations. It is worth mentioning that the same equa- 
tions are obtained by performing a rigorous derivation 
based on the multiscalc expansion technique, which we 
shall publish elsewhere. Notice that Eqs. H]) have only 
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two independent parameters: pump A and quotient of 
cavity losses a, which we take a > 1 without loss of gen- 
erality. 

Eqs. ([4]) have four possible stationary ho- 
mogeneous solutions: The laser-off solution, 
F\ = i*2 = 0, the strong unidirectional (or single- 
mode) solution {i 7 ! = V 'A — 1, F^ — 0}, the u>eafc 
unidirectional solution {i<\ = 0,F 2 = -\A4 — <r}, 
and the bidirectional (or bimode) solution 
{Fx = (A — 2a+l) /3, F 2 = (A - 2 + a) /3}. The 
linear stability analyses of these solutions shows that: 
(i), The laser-off solution losses its stability at A = 1 
giving rise to the strong unidirectional solution, which 
is always linearly stable; (ii), the weak unidirectional 
solution is linearly stable for A > 2a — 1; and (iii), 
the bidirectional solution is always unstable. Thus 
bidirectional homogeneous cw emission is impossible, as 
it is well known from the plane-wave theory Q. Our 
basic hypothesis is that stable bidirectional cw emission 
could be possible in a nonhomogeneous solution and this 
is the possibility we examine now through the numerical 
integration of Eqs. ([4]) with the purpose of finding CS 
solutions. 

We shall concentrate here in the one-transverse di- 
mension limit, i.e., we assume that diffraction along one 
transverse direction is avoided by some means, while it 
acts freely on the other (orthogonal) coordinate (think, 
e.g., of a waveguide configuration in one transverse di- 
rection) . We have used a split-step algorithm for solving 
Eqs. ([U with a grid of 512 spatial points and a total 
integration length L/y/a = 1 with a = 1CP 5 . In order to 
approach realistic conditions, here we consider that pump 
is limited in the transverse direction and further assume a 
top-hat like profile for A (for computational convenience 



with £ p = 0.45 and 



we take A(£) = Aexp - (£/£p) 2 " 

n = 10). Nevertheless, we have also carried out integra- 
tions with an unlimited pump function in order to check 
the robustness of our results. 

We have first investigated the dynamics of the system 
when a real positive pulse is injected in the weak mode 
(its maximum value equals the weak-mode homogeneous 
solution), and the strong-mode homogeneous solution, 
with a dip arriving to zero in its center, is injected in the 
strong mode. Remarkably, the results are independent 
of the exact shape and width of the pulse (which are the 
same as that of the dip). 

For o~ = l this initial condition relaxes to a final state in 
which the weak mode is null in all the space irrespective 
of the pump value [9(. A quite different result is reached 
as a is made different from unity and A > 1.05: a narrow 
domain of emission in the weak mode, and a correspond- 
ing dip in the emission of the strong mode, are formed 
where the initial condition was injected, giving rise to a 
localized structure. 

In Fig. 1 we show, in the (a, A) plane, the domain 
of existence of these localized structures. In the inset 
an example of the intensity profiles of the two fields is 



shown. The bright CS of the weak mode is nearly Gaus- 
sian, while the dark CS -or, more precisely, the grey CS, 
as the minimum intensity is clearly different from zero- 
of the strong mode corresponds to a nearly supergaus- 
sian dip on the strong-mode homogeneous solution. The 
CSs represented in the inset of Fig. 1 correspond to the 
domain in which they are stationary but the existence 
of a Hopf bifurcation as a is decreased must be noticed. 
The Hopf bifurcation is supercritical, i.e. when reached 
(decreasing a) small oscillations appear that affect both 
the height and width of the CS, oscillations that grow in 
amplitude as a is further decreased. In Fig. 2 an example 
of the dynamic behavior exhibited by the CS is shown. 
Fig. 2(a) shows the spatiotemporal oscillation of the two 
fields, the maximum of weak field, F2, exhibits periodic 
oscillations, see Fig. 2(b). 

The continuous line in Fig. 1 corresponds to the 
boundary where the CSs cease to exist. When cross- 
ing this line, the bright CS of the weak mode simply 
decreases in amplitude, quickly, till it disappears, giving 
rise to homogeneous emission in the strong mode. Some- 
thing similar occurs at a — 1, where the oscillating CSs 
collapses. 

After establishing the existence of these CSs, an im- 
portant issue is to determine whether they can be ma- 
nipulated, i.e., how can they be written/erased. Notice 
that in the bidirectional laser one can inject signals in 
either the strong or the weak mode or in both. The re- 
sults we have just commented are obtained from an ini- 
tial condition which, to some extent, can be understood 
as the simultaneous injection of a positive pulse in the 
weak mode and a negative pulse in the strong one. But 
from an experimental perspective this could be difficult 
to implement. Moreover, one must take into account the 
unavoidable presence of noise in real devices. 

We have thus carried out a series of numerical integra- 
tions in which we have tested the writing and erasing of 
CSs by injecting appropriate (i.e., with the right phase) 
pulses in only one of the two counterpropagating fields 
in the presence of additive noise (random noise with am- 
plitude 10~ 5 ) in Eqs. (U]). The results are summarized 
in Fig. 3. The initial condition consists of emission only 
in the strong mode. We start with homogeneous steady 
state emission in the strong mode. At t = 0, a positive 
Gaussian pulse is added only in the weak mode for a short 
time, the result being the formation of the CSs. The CSs 
are then removed by injecting a positive Gaussian in the 
strong mode, at r = 2-10 4 . Then, at r = 2.5T0 4 a neg- 
ative Gaussian is injected in the strong mode, and the 
CSs are formed again, which are later removed by inject- 
ing again at r = 3. 5 TO 4 a Gaussian pulse in the weak 
field but with the opposite phase (the injection times are 
marked in the figure with a horizontal white line) . These 
numerical results prove the robustness and addresability 
of these localized structures. 

We have also assessed the possibility that the system 
supports clusters formed by an array of CSs. The numer- 
ical results ensured that several CSs can be switched on 
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and that this solution is stable and robust in the pres- 
ence of noise. Figure 4 shows the spatiotemporal evo- 
lution where a cluster formed by four CSs is developed 
(a) from the strong mode emission initial condition when 
four Gaussian pulses are injected in the weak field and 
the final intensity profile (b). 

In conclusion we have shown that cavity solitons form 
in class-A bidirectional lasers when the two counterprop- 
agating fields experience slightly different cavity losses. 
Wc find that it is particularly remarkable the fact that (i) 



the CSs in the two fields are complementary, i.e., when 
one field exhibits a bright CS, the counterpropagating 
field shows a dark CS in the same position, and (ii) the 
CSs can be written/erased by acting on only one of the 
two fields. 
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FIGURE CAPTIONS 

Fig.l. Bifurcation diagram of CSs in the < a, A > 
plane. Continuous and dashed lines represent the bound- 
aries of existence of stationary and dynamic CSs, respec- 
tively. The inset shows a stationary CS obtained for 
a = 1.01 and A = 1.4. Continuous (dashed) line cor- 
responds to the field amplitudes in the weak (strong) 
modes. 

Fig. 2. Pulsing CS obtained for a = 1.005 and A = 
1.4. The spatiotemporal evolution of the pulsing CS is 
shown in (a), where r runs from 2.9 10 5 to 3 10 5 and £ 
from —0.5 to 0.5, the maximum of the CS in the weak 
field oscillates periodically in time (b). 

Fig. 3. Switching on and off of CSs by different tech- 
niques. Parameters are a = 1.015 and A = 1.4, £ running 
from —0.5 to 0.5. The white lines indicate the times at 
which pulses are injected. For details, see text. 

Fig. 4. Cluster formed by four CSs for a = 1.005 and 
A = 1.4. Time r runs from to 2.5 10 4 . Spatiotemporal 
evolution (a) from the injection and final intensity profile 
(b). 
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